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Robust  C(a)-type  Tests 
For  Linear  Models 


by 

Elvezio  Ronchetti 
Department  of  Statistics 
Princeton  University 


ABSTRACT 


A  new  class  of  tests  that  can  be  viewed  as  a  generalization  of  Neyman's 
optimal  C(a)  test  is  introduced.  An  optimally  robust  test  which  maximizes 
the  asymptotic  power  within  this  class,  subject  to  a  bounded  influence 
function.  Is  selected.  It  Is  shown  that  it  is  equivalent  to  a  certain 
asymptotically  minimax  test  proposed  by  Wang  (1981).  Finally,  some  numerical 
results  on  the  asymptotic  behaviour  of  robust  C(a)-type  tests  under  several 
errors'  and  carriers'  distributions  are  discussed. 


AMS  (1980)  SUBJECT  CLASSIFICATION:  Primary  62F35 

Secondary  62F03. 


KEY  WORDS  AND  PHRASES:  robust  testing.  Influence  function,  asymptotic  power, 

C(a)  tests,  optimally  robust  tests,  regression  model. 


1.  INTRODUCTION 


Let  {(x^.y^)  :  1*1  .....  n>  be  a  sequence  of  Independent  Identical 
distributed  random  variables  such  that 

*  x]e  +  e^  ,  1*1  .....  n  ,  (1.1) 

where  y^  is  the  1th  observation,  x^eFp  Is  the  1th  row  (written  as  a 
column  vector)  of  the  design  matrix,  6  e  ©C-Fp  a  p-vector  of  unknown 
parameters  and  e^eF  the  1th  error.  Suppose  that  e^  Is  Independent  of 
x.j  and  Is  distributed  according  to  a  normal  N(0,a  )  .  Moreover,  denote 
by  K(x)  the  distribution  of  the  x's  with  respect  to  some  o-flnite 
measure  y  and  by  FQ(x,y)  the  joint  distribution  of  (x^y^)  . 

In  the  classical  parametric  approach,  there  are  several  procedures 
one  can  use  for  subhypothesis  testing:  the  F-test,  which  is  equivalent 
to  the  likelihood  ratio  test,  the  Wald  test  and  the  C(a)  test. 

These  tests  are  asymptotically  equivalent  and,  when  the  errors  are  normally 
distributed,  they  are  optimal;  see,  for  Instance,  Cox  and  Hinkley  (1974). 
However,  these  test  procedures  suffer  similar  robustness  problems  as  the 
least  squares  estimators.  Although  they  are  moderately  robust  with  respect 
to  the  level,  they  do  lose  power  rapidly  in  the  presence  of  small  departures 
from  the  normality  assumption  on  the  errors;  cf.  Hampel  (1973,  1978),  Schrader 
and  Hettmansperger  (1980).  Robust  alternatives  to  the  F-test  and  the  Wald 
test  have  been  Investigated  recently  by  Schrader  and  McKean  (1977),  Schrader 
and  Hettmansperger  (1980)  and  Ronchetti  (1982a,  1982b,  1984). 
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In  this  paper  we  focus  on  the  class  of  C(ct)  tests.  In  section  2 
we  define  a  new  class  of  tests  which  generalizes  the  optimal  C(a)  tests 
Introduced  by  Neyman  (1958,  1979).  In  section  3  we  investigate  their 
robustness  properties  by  means  of  the  Influence  function  and  we  compute 
their  asymptotic  power.  This  allows  us  to  select  the  optimally  robust 
C(o)-type  test,  that  Is  a  test  which  maximizes  the  asymptotic  power 
within  this  new  class  (efficiency  condition),  subject  to  a  bounded 
influence  function  of  the  test  statistic  (robustness  requirement).  In 
section  4  we  show  the  equivalence  between  our  optimally  robust  C(a)-type 
test  and  an  asymptotically  minimax  test  Introduced  by  Wang  (1981).  This 
points  out  the  strong  relation  between  the  minimax  approach  and  our 
approach  to  robust  testing.  Finally,  section  5  presents  some  numerical 
results  on  the  asymptotic  behaviour  of  C(a)-type  tests  under  several 
distributions  of  the  errors  and  of  the  carriers. 


Consider  the  regression  model  (1.1).  Though  C(a)  tests  can  be  defined 
for  testing  more  than  one  linear  hypothesis  on  e  (see  BUhler  and  Purl, 
1966),  we  focus  here  on  the  generalization  of  Neyman’s  original  definition 
(see  Neyman,  1958),  that  Is  we  shall  Introduce  a  new  class  of  procedures 
for  testing 


where  0^  denotes  the  jth  component  of  the  vector  0  .  For  a  given 
p-vector  x  ,  we  denote  by  Its  "nuisance  part",  that  Is 

X(l)  *  (x<l> . x^V  and  by  x^j  Its  last  component  x^  . 

Moreover,  let  x  *  •  We  denote  by  ,  M^)  *  M(21)  '  M(22) 

the  submatrices  of  a  (p-l)xl  partition  of  a  pxp  matrix  M  .  Note 
that  Mj22)  Is  simply  mpp  ,  the  (p,p)  component  of  M  . 

The  class  of  tests  we  shall  define  Is  based  on  a  function 

n:FP  xF+F  ,  (x.r)  -*■  n(x,r)  , 

which  satisfies  the  following  conditions: 

(2.ETA1)  (1)  n(x,0  Is  continuous  and  odd  for  all  x  e  , 

(11)  ti(x.r)  >0  for  all  x  e  F  ,  r  e  F+ 

(2.ETA2)  n(x,*)  Is  differentiable  on  R\D(x;n)  for  all  x  c  , 
where  P(x;n)  Is  a  finite  set. 


Let  n'(x.r)  :=  (3/3  r)n(x,r)  If  xeRp  ,  reR\0(x;n) 

:*  0  otherwise. 

and  assume: 

supr|n'(x,r)|  <•  for  all  xeFp  . 

We  shall  also  assume  the  following  regularity  condition: 

(2.ETA3)  M  :«  En'(x,r)xxT  exists  and  Is  nonsingular. 

Definition  2.1  The  class  of  C(a)-type  tests  for  linear  models  Is  defined 
by  means  of  the  test  statistics 

Z  (6;n)  :*  n-**  I  Ti(xi,(yi-xl9)/o)-(U"1xi)^p^ 

n  1-1  1  1  1  1  (2.2) 

■  ""\p  Ji  n(xi>(yrxf5)/0)'[:>iP)-u(2i)'U(n)(xi>(i)] 

where  n  satisfies  the  conditions  (2.  ETA1),  (2.ETA2)  and  U  Is  the  lower 
triangular  matrix  with  positive  diagonal  elements  such  that  UUT  -  M  .  "Large" 
values  of  Zn  are  significant. 

Remark  1.  C(a)  tests  were  Introduced  by  Neyman  (1958)  for  a  general  parametric 
model  and  extended  to  the  robust  testing  problem  by  Wang  (1981,  Remarks  2  and 
3,  p.  1100)  who  was  able  to  derive  an  (asymptotically  minimax)  robust  version 
of  the  optimal  C(a)  test.  We  shall  discuss  the  connections  with  this  test 
In  section  4.  If  we  put  n(x,r)  -  -♦,(r)/4(r)  *  r  in  (2.2),  where  <t> 
is  the  standard  normal  density  (the  density  of  the  error's  distribution),  Zn 
becomes  the  test  statistic  of  the  optimal  C(a)  test  obtained  by  Neyman  (1958, 
p.  228).  In  this  sense  the  tests  defined  by  (2.2)  can  be  called  C(a)-type  tests. 


Remark  2.  The  test  statistic  Zn  depends  on  the  unknown  nuisance  parameters 
and  on  the  scale  parameter  o  .  In  section  3  we  shall 
discuss  the  properties  of  a  studentlzed  version  of  Zn  .  From  now  on  we 
assume  for  simplicity  o  *  1  . 

Remark  3.  In  order  to  determine  asymptotic  critical  regions  of  the  test, 
we  shall  compute  the  asymptotic  distribution  of  ZR  In  section  3. 
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3.  OPTIMALLY  ROBUST  C(ct)-TYPE  TESTS 

In  this  section  we  investigate  the  robustness  properties  of  C(o)-type 
tests  by  means  of  the  influence  function  and  we  select  the  optimally  robust 
test  In  this  class.  This  Is  a  test  that  maximizes  the  asymptotic  power  within 
the  class  of  C(a)-type  test,  subject  to  a  bounded  Influence  function. 

The  notion  of  Influence  function  for  estimators  was  introduced  by 
Hampel  (1968,  1974).  It  is  essentially  the  first  derivative  of  an  estimator 
viewed  as  functional,  and  describes  the  normalized  Influence  of  an  infinitesimal 
observation  on  the  estimator.  Formally,  suppose  the  estimator  Tn  can  be  ' 

written  as  functional  T  of  the  empirical  distribution  function  , 

Tn  *  •  Then  the  Influence  function  of  T  at  F  is  given  by 

t 

IF(z;T,F)  «  HiV^(T((l-e)F+e62)  -  T(F))/e  ,  (3.1) 

where  6Z  is  the  distribution  that  puts  mass  1  at  the  point  z  .  It  is  ! 

then  clear  that,  from  a  robustness  point  of  view,  boundedness  is  a  desirable 
property  of  IF. 

The  concept  of  Influence  function  has  been  extended  to  tests  by  Rousseeuw  ! 

and  the  author,  and  independently  by  Lambert;  see  Ronchetti  (1979,  1982a); 

Rousseeuw  and  Ronchetti  (1979,  1981);  Hampel,  Ronchetti,  Rousseeuw,  Stahel 
(1984);  Lambert  (1981).  It  turns  out  that  the  influence  function  defined  on  ' 

the  test  statistic  is  proportional  to  the  Influence  of  an  infinitesimal  observa¬ 
tion  on  the  level  and  on  the  power  (or,  as  in  Lamber  (1981),  on  the  P-value) 
of  the  test.  Therefore,  a  test  statistic  with  a  bounded  influence  function  i 

guarantees  robustness  of  validity  and  robustness  of  efficiency  for  the  test. 

< 
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Let  us  now  compute  the  influence  function  of  our  C(a)-type  tests.  Let  Z 
be  the  functional  corresponding  to  the  test  statistic  Zn  ,  that  is 

Z(F)  */n(x,y-xT0)*(U"1x)^p^dF(x,y)  . 

Then  Z(F^)  *  n_lsZn  ,  where  F^  is  the  empirical  distribution  function. 

Proposition  3.1  Let  Fg  be  the  model  distribution  under  the  null  hypothesis. 
If  (2.ETA1)  and  (2.ETA2)  hold,  the  influence  function  of  Z  is  given  by 

IF(x,y;Z,Fg)  =  n(x,y-xT6)  •(lf1x)^  (3.2) 

Proof.  The  result  follows  easily  applying  the  definition  (3.1)  to  Z  and 
noting  that  by  (2.ETA1)  (i)  En(x,r)x  *  0  . 

From  (3.2)  we  see  that  Heyman's  optimal  C(a)  test  (n(x,r)  *  r)  ,  though 
asymptotically  efficient,  has  an  unbounded  influence  function.  Our  goal  will 
be  to  compromise  between  efficiency  and  robustness  within  the  class  of 
C(o)-type  tests  by  finding  a  test  that  maximizes  the  asymptotic  power,  under 
a  bound  on  the  Influence  function.  Let  us  now  compute  the  asymptotic  power 
of  a  C (a) -type  test  defined  by  rj  . 

Proposition  3.2  Besides  (2.ETA1),  (2.ETA2)  and  (2.ETA3)  assume 


(3.AS3)  En^(x,r)  xj^  <  ® 


» 


where  |*|  denotes  the  Euclidean  norm.  Put  Xp  :«  En2(x,r)*|(U”1x)^|2  . 
Then,  under  the  sequence  of  alternatives 

H(n)  :  0  *  0  +  n->sA  ,  (3.3) 

where  A  *  (0,...,0,A^)T  , 

the  test  statistic  Zn  has  asymptotically  a  normal  distribution  with  mean 
Upp-A^  and  variance  Xp  .  Moreover,  the  asymptotic  power  of  the 
C(a)-type  test  at  the  level  o  defined  by  Zn  is  given  by 

1  -  *(,‘1»-o) -A(p)-uppA^)  , 

where  *  Is  the  cumulative  standard  normal. 

Proof.  It  suffices  to  compute  E„Z_  under  the  sequence  (3.3).  Define: 

■  v  n 

E(0)  :«  /n(x,y-xTe)-(U'1x)^p^f0(x,y)-dy(x)dy  , 

where  fQ(x,y)  *  <f(y-xTe),k(x)  .  Then,  for  j  ■  l,...,p-l  ,  we  have 

C(3/30^)S(0)3g  *  -/n'(x,y-xT0)x^'(lf  M^dFgfx.y) 

+  /n(x,y-xTe)-(U”1x)^p^C(3/30^h^e(x,y)3gdy(x)dy  . 


(3.4) 
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Using 


ar'xi'p’  •  <'pJ-(-u(21)U(}1)x(1)+x(p>)  , 


(3.5) 


we  get 


En'  (x,r)x^1}  •  (u"3*)  (p) 


(3.6) 


upp"  (”U(21)  *u(ii)  *Eti '  ^x,r) x(1) +  EnMx^x^x^) 
upp‘<"u(2i)%uai),uai)‘ua1)  +Ma2))  *  o  , 

and  the  first  term  of  the  right  hand  side  of  (3.4)  vanishes.  In  a  similar 
way  one  can  prove  that  also  the  second  term  of  the  right  hand  side  of  (3.4) 
equals  0.  Moreover, 


[(3/36 *p) )  U6)  ]g  - 

Jn(x,y-xT6)  •  (U_1x)  (p)  •  I(3/3e{p))fe(xry)]^dy(x)dy 
En (x,r)  •  (u”*x)  ^  *r*x^p^  * 

and  by  Integrating  by  parts, 

En'(x,r)*(«‘1x)(p).x(p)  - 


UPP‘<‘D(21)*U'(U)*M(12)  +  "pp>  *  °pj«“pp>  *  “pp 


(3.7) 
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Finally,  denoting  by  a  dot  the  differentiation  with  respect  to  0  ,  we  obtain 

<$<6))(j>«0  ,  for  j  (3.8) 

<l(8))(p)  «  Upp  (3.9) 

and  by  (2.ETA1)  (1) 

E(0)  =  0  . 

Therefore,  using  a  Taylor  expansion  we  get 

E0zn  -  n1/2$<0)  -  n1/2ST<e)*(e~e)  +  ©(||e-e||) 

*  uppA*p)  +  o(n-1^2)  , 

This  completes  the  proof. 

Remark  1.  The  test  defined  by  the  test  statistic  ZR  depends  on  the  unknown 

nuisance  parameter  0  .  By  techniques  similar  to  those  used  by  Wang  (1981) 

one  can  show  that  the  result  of  Proposition  3.2  holds  If  we  substitute  0 
k  ' 

by  a  suitable  (n  -consistent)  estimate  Tn  (see  Wang,  1981,  p.  1099). 
Moreover,  the  result  of  Proposition  3.1  still  holds.  To  see  this,  suppose 
that  the  influence  function  of  T  exists  and  define 

2(F)  -  /n(x,y-xT-T(F))-(U_1x)(p)dF(x,y)  . 


(3.10) 


IF(x,yjZ,Fg)  -  IO/3e)Z((l-e)Fe  +  £6(x,y)>le-0 

■  n(x,y-xT0)  •  (U-1x)  ^  +  Jtj(«#v-sT8)  •  (U^s)  (p*dFg(s,v) 
-c/n* (8,V-ST6) • (u”1*) *P* •»TdF§(s,v) ) *IF(x,y*T  ,  Fj)  . 

Now  the  second  term  of  the  right  hand  side  of  (3.10)  Is  equal  0  In  view 
of  (3.5)  and  (2.ETA1)  (1),  and  so  Is  the  third  term  by  (3.6). 

The  next  Proposition  gives  the  optimally  robust  C(a)-type  test. 

Proposition  3.3  Under  (2.ETA1),  (2.ETA2),  (2.ETA3),  (3.AS3),  the  test  which 
maximizes  the  asymptotic  power  within  the  class  of  C(a)-type  tests,  subject 
to  a  bounded  Influence  function,  is  defined  by 

n0(x,r)  =  |z(p)r\;(r|z(p)|)  *  i|>  (p)  (r)  • 

c/|z'p'| 

where  z  *  U^x  ,  UQ  is  the  lower  triangular  matrix  defined  implicitly  by 
the  matrix  equation 

E(2*(c/|z(p)|)-1)zzT  -  I  , 

♦c(t)  *  m1n(c,max(t,-c))  Is  the  Huber  ^-function,  and  c  is  a  given  positive 
constant  depending  on  the  bound  on  the  Influence  function. 


2 

According  to  Proposition  3.1  and  3.2  we  have  to  maximize  Upp/Xp  »  subject 
to  a  bound  on  the  right  hand  side  of  (3.2).  The  same  problem  has  been  solved 
In  Ronchettl  (1984;  Theorem  4.3)  In  connection  with  the  class  of  so-called 
x-tests.  Therefore  we  refer  to  that  paper  for  the  complete  proof;  cf.  also 
Ronchettl  (1982a,  Theorem  3). 
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4.  CONNECTION  WITH  AN  ASYMPTOTICALLY  MINIMAX  TEST 


In  this  section  we  describe  the  connections  between  the  optimally  robust 
C(a)-type  test  and  a  test  introduced  by  Wang  (1981).  Wang  studies  the  testing 
problem  using  minimax  techniques.  He  considers  the  following  situation. 

Suppose  we  are  given  a  parametric  model  CFQ  :6  e6  C  Fp}  and  suppose 
we  want  to  test  a  hypothesis  on  one  component  of  the  parameter  e  ,  say 
6^  *  0  ,  the  other  (p-1)  components  being  nuisance  parameters.  Then, 
using  the  technique  of  shrinking  neighbourhoods,  Wang  Is  able  to  derive  an 
asymptotically  minimax  test  in  the  case,  where  the  model  distribution  is 
Indexed  by  nuisance  parameters.  Let  us  write  this  result  in  the  situation 
of  the  linear  model.  Define: 


r  *« 


y  -  x  e 


a  :• 


u(1> . .«*>>*  . 


d(x,a(ij)  i-  x*p)  *  x<P>  +2j”J  a^x*^  , 

A(atl)  ,6)  «  A*(x,y  ;a(1j  ,6)  :»  r*d(x,a )  .  (4.1) 

For  a  given  e  >  0  ,  6^  >  0  ,  define  vo*a(i)'®^  1111(3  Vl*a(l)'®^ 

implicitly  by  means  of  the  equations 

E(max{(A(a(1)  ,6)  -  Va(1)  ,6) )  ,0})  -  e/6x  (4.2) 

Elmax{(V0(Sjlj  ,6)  -  A(a^^  ,0))  ,0} +A(atl^  ,6)1  ■  c/6 A  .(4.3) 

Moreover,  let  a^  ■  (6)  be  the  solution  to  the  equation 

_  V.  (am  ,0) 

•  r.*(l))  -0  ,  (4.4) 


and  define 
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o  *  «.V«(cX/.x»0) 

v2(8)  :*E(t*(a(l).e)JvJ(o«)t5),2  (4.5) 

Proposition  4.1  (Wang,  1981,  p.  1099,  p.  1104)  The  test  defined  by  the 
critical  region 

(Yn(0)  >  r^l-o)  +  eV1(o(1)(8),e)/v(0)}  ,  (4.6) 


where 

Yn(e)  :*  n  Mv(e))  CA*(xi,y1;o(1)(e),e)DvJ(aji1j^)  , 

Is  an  asymptotically  minimax  test  at  level  a  . 

The  test  defined  by  (4.6)  will  be  called  Wang  test.  In  order  to  have 
a  better  performance  at  the  model,  Wang  (1981,  p.  1105)  proposes  a  modification 
of  the  test  (4.6)  and  defines  a  test  by  means  of  the  following  critical  region 

(Yn(e)>*"1(l-a)}  (4.7) 

The  test  defined  by  (4.7)  will  be  called  modified  Wang  test.  Then  we  have 
the  following  result. 

Proposition  4.2  The  modified  Wang  test  Is  a  C(a)-type  test.  It  is  equivalent 
to  the  optimally  robust  C(a)-type  test  defined  by  ns  (given  by  Proposition  3.3). 
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^  •* 

Proof.  We  prove  that  Yn(e)  *  (n0) •zn(e;n0)  •  First,  consider 

equation  (4.4).  Using  (4.1)  we  have 

vo  v, 

0  *  EnA]V0*r*x(l)>  “  EUr«dJv1*r«x(1)) 

f  V./|d| 

-  /  <Kx,a{1))«x(1)-(/[r]vJ/|dj.rd4(r))dK(x)  , 

and  (4.4)  becomes 

/(♦<VI'au.«,x)>  |)-*(v#/|4(x,«a)>  l>) -acx.ajDlJCd)"1**  "° 

(4.8) 


Now,  combining  (4.2)  and  (4.3)  and  noting  that  EA  -  0  ,  we  obtain 

/(A-v1).1Ujv1,  ♦  /<*-V0'-l(A*v0)  -  0  * 

and  performing  these  integrations  we  get  finally  by  (4.8) 

VQ  *  -Vx  («:  -V)  .  (4.9) 

For  a  given  positive  constant  c  ,  define  ng  and  UQ  as  in  Proposition  3.3. 
Moreover,  choose 


» ■  c-eVpp  • 


(4.10) 


16 


Then 

“(1)  *  “*V(21)’^U0*(llp  1  *4-U) 

and 

d(x,a(i))  s  x^"(uo)(2if  ((ucP(ii))  lx{i) 

(4.12) 

=  (U51x><P>.(l^)pp  • 

(The  left  member  of  (4.8)  becomes 

EC(2*(c/|(U"1x)(p) |)-l)*(U"1x)(p)-(U0)pp-x(1)D 
*  tEn^(x,r)-(U^)(p)-x{1)3-(U0)pp  -  0  by  (3.6)  .) 

Moreover,  by  (4.5),  (4.9),  (4.10)  and  (4.12)  we  have 

vJ(6)  -  E(tr.dl^)2  .  B(aJ(*i«(l))-*J/|(S|(r)) 

-  (D0,PP,E<^/|J (p)  ,r)  •  l!<P>  l2>  (4.13) 

■  ‘Vpp-yno* 


9 
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where 


»  (U^x)^  and  Ap(n0)  ■  En^x.r)- |z^p^ |2  . 


Finally  we  get 


ve» 


Si  .  n-1'2. 


((0o,pp-xp/j<no,rl-1wlr-d,-v 


n-l/2.x-l/2(V  •ti;i((U0)pp)-1.*e(Vpp(r.*<P)  •  (U0)pp) 
n-V2  .  ^-1/2  (t,q  , .  Z^c  (r •  z (p> ) 


Xp1/2(VV5’V  • 


This  completes  the  proof. 


Proposition  4.3  The  Wang  test  is  always  less  efficient  at  the  model  than 
the  modified  Wang  test. 

Proof.  The  modified  Wang  procedure  has  the  same  asymptotic  power  as  the 
optimal  robust  C(a)-type  test  defined  by  n0  •  Therefore,  the  square  of 
Its  efficacy  equals  (U0)ppAp(n0)  . 
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On  the  other  hand,  the  square  of  the  efficacy  of  the  Wang  test  Is  given 
by  the  formula  (see  Wang,  1981,  p.  1104) 


s2(e)  -  (v(e)  -  (c^W-vj/vfe))2 


a 

d 


Thus,  the  relative  efficiency  of  the  Wang  test  with  respect  to  the  modified 
Wang  test  can  be  computed  as 


eff{Wang  test,  modified  Wang  test} 


s2<e)/((U0)2pAp(n0)> 


hence,  using  (4.10)  and  (4.13) 


[1  (nJ/lV, J  * 

P  0  0  PP 


•  xl/2  <no} + <e/5i)  * c*  <V  pp/  ( (V  pp* Xp  (n())  * 1 


WPP'Ap 


(X  (no)  ♦  (e/6x)-c/(u0>pp)  • 


Using  (4.2)  we  get  finally 


eff{Wang  test,  modified  Want  test} 

*  l-c2EC(|z{p)|/c)-*'(c/lz(p)|)  -*(-c/|z(p)|)3  <  1 


This  completes  the  proof. 
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5.  ASYMPTOTIC  BEHAVIOUR 


In  this  section  we  present  some  numerical  results  on  the  asymptotic 
behaviour  of  C(a)-type  tests  under  several  distributions.  In  the  exposi¬ 
tion  we  apply  the  methods  which  have  been  used  by  Maronna,  Bustos  and  Yohal 
(1979)  for  comparing  different  regression  estimators.  We  consider  simple 
regressslon: 


y  *  x  6  +  e  , 

where  x  -  (l,x^)T  ,  6  «  {6^  ,  e^)T 
We  want  to  test  the  hypothesis 


Hq  :  «  0  . 


Let  G  and  K  be  the  distributions  of  e  and  x^  ,  respectively. 

For  x<«  and  e  we  choose  the  following  contaminated  normal  distributions 


(l-e)*(*)  +  e*(7  s) 


with  the  following  parameters 


e  :  0  0.1  0.1  0.05 

s  :  1  3  5  10 


We  want  to  compare  the  asymptotic  behaviour  of  the  optimally  robust  C(a)-type 
test  (see  Proposition  3.3) 

n0(x.r)  -  ♦Ci/|2(2)|(r) 

with  that  of  another  C(a)-type  test  with  bounded  Influence  function,  namely 

n^x.r)  *  wc  ( Jz^2^ | )ipc  (r)  ,  (5.1) 

2  3 

where  Wc(|t|)  *  \pc(t)/t  *  min(l ,c/ j tj )  .  The  C(a)-type  test  defined 
by  (5.1)  Is  the  solution  of  the  optimality  problem  considered  In  Proposition 
3.3  within  the  restricted  class  of  functions  n  satisfying  n(x,r)  «  w(x)-ip(r) 
see  Ronchetti  (1982b).  The  subscript  "M"  reminds  that  is  a  function  of 
Mallows'  form. 

In  order  to  Investigate  the  behaviour  of  a  C(a)-type  test  defined  by 
a  "redescending"  n-functlon,  we  consider  also  the  following  procedure: 

Tin|j(x,r)  *wc^(|z^  ^ | )*(|'(r;Cg,ic*,Cg,Ag,Bg)  , 

where  is  defined  by 

(i(r»Cj»ic*iCg,ftgiBj)  ■  r  If  0  ^  | r |  ^  Cg 

■  o-tanhtB(Cg-|r|)Ds1gn(r)  cg<  |r|  <  cg 

■  0  c5<  |  r  J  , 
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o  *  (Ag(ic*-1 ) )**  ,  g  «  ((<*-l)Bg/Ag)Js/2  .  The  function  $  defines  the 
socalled  "hyperbolic  tangent  estimator"  for  location;  see  Hampel.  Rousseeuw, 

Ronchettl  (1981).  Note  that  Cg.Ag.Bg  are  computed  Implicitly  In  terms  of 
c5  and  k*  . 

For  each  test  we  compute  the  standardized  sensitivities  at  the  normal 

model 

u22*supx,rlr,(x *r>  I  ‘  lz^  I 

and  the  efficacies 

2 

U22^X2 

under  several  distributions. 

Note  that  the  efficacy  of  the  tests  defined  by  and  factorizes; 
in  these  cases  we  have 

u22^2  *  DX‘DR  * 

2  2 

where  DX  depends  only  on  wc  and  K  and  DR  =6  /A  .with  A  *  E\pc  .  B  *  . 

The  constants  are  chosen  so  that  all  the  tests  have  the  same  asymptotic 
efficiency,  when  ~  N(0,1)  and  e  ~  N(0,1)  .  ' 

Table  1  describes  the  calibration  as  well  as  the  standardized  sensitivities 
of  the  tests.  From  Table  2  one  can  obtain  the  following  conclusions: 

( 


4 
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11  "HT  Is  better  than  for  all  distributions  under  consideration; 

2)  t)q  Is  better  than  when  the  distribution  of  the  errors  has 

moderate  tails; 

3)  tiq  has  the  better  standardized  sensitivities  (computed  at  the  normal 
model ) . 

Table  3  gives  the  factors  for  the  and  r^-test. 
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Table  li  Calibration  and  standardized  sensitivities 
at  x(2)  ~JT{0,1>  ,  e  ~*(0,1) 


Table  2:  Asymptotic  efficacies 
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K 

e-0 

Ml 


M  *1 
8-3 


e-  *1 
MS 


e*  *05 

s-lo 


TEST 


e-o 

••l 

e»  .1 
•■3 

€■  .1 
8-5 

.95 

.75 

.71 

.95 

.77 

.77 

.95 

.78 

.72 

1.52 

1.21 

1.13 

1.53 

1.23 

1.23 

1.54 

1.27 

1.16 

2.44 

1.94 

1.81 

2.44 

1.97 

1.97 

2.51 

2.06 

1.86 

2.99 

2.38 

2.23 

3.00 

2.42 

2.42 

3;  24 

2;bi 

2;  32 

C-  .05 

ST. SENS 

8-10 

•  80 

3.16 

.87 

3.37 

•7B 

2.87 

1.28 

2.70 

1.40 

2.91 

1.25 

2.34 

2.04 

2.33 

2.24  . 

2.51 

2.00 

1.94 

2.51 

2.23 

2.76 

2.38 

2i50 

1.85 

Table  3*  Factorization  for  nMT-teat 


C"0 

■■1 

t*  .1 
•■3 

e-  .1 
a-5 

€•  .05 
••10 

DR 

•970 

.771 

.722 

.614 

•973 

•765 

.785 

.694 

OX 

.979 

1.569 

2.512 

3.067 

